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ABSTRACT 



We present a measurement of the quasar luminosity function in the range 0.68 < z < 4 down to extinction corrected 
magnitude gdered = 22.5, using a simple and well understood target selection technique based on the time-variability 
of quasars. The completeness of our sample was derived directly from a control sample of quasars, without requiring 
complex simulations of quasar light-curves or colors. A total of 1877 quasar spectra were obtained from dedicated 
programs on the Sloan telescope (as part of the SDSS-III/BOSS survey) and on the Multiple Mirror Telescope. They 

■ allowed us to derive the quasar luminosity function. It agrees well with previously published results from Croom et al. 

■ (2009) in the common redshift range 0.68 < z < 2.6. Our deeper data also allow us to extend the measurement to z — 4. 
£C} ' We measured quasar densities to gderod < 22.5, obtaining 30 QSO per deg 2 at z < 1, 99 QSO per deg 2 for 1 < z < 2.15, 

and 47 QSO per deg 2 at z > 2.15. Using pure luminosity evolution models, we fitted our LF measurements combined 
with the data from Croom et al. (2009), and predicted quasar number counts as a function of redshift and observed 
magnitude. These predictions are useful inputs for future cosmology surveys such as those relying on the observation 
of quasars to measure baryon acoustic oscillations. 

Key words. Quasars: general, dark energy, surveys 



X. 

H . 1- Introduction 



Riehards et all. 12001 12004 120091: ICroom et all 120091: 
Bow et al. 1. 120111 120121) present serious drawbacks for the 



The measure ment of^ the B a ryon Acou stic Oscillation selection of quasars at redshifts near z~2.7, which oc- 

(BAO) scale flEisenstein et al.|, |2005|, |2007|) relies on large cupy s i m il ar regions of ugriz color-space as the more 

samples of objects selected with an unbiased method. To numerous white dwarfs and blue halo sta r s (iFanl . 119991: 

probe the distant Universe, quasars appear to be one of the |Rj c hards et~aIL l200l IWorseck fc Prochaskd I2010D. To cir- 



sources of choice, since they are both among the brightest cumvent this dimcultvT lPalanaue-Delabrouille et all (|201ll) 



extragalactic objects, and expected to be present at sum- developed a selection algorithm relying on the time variabil- 

ciently high density. ity of quasar fluxes. This technique was te sted in 2010 as 

The selection of quasars to redshift z ~ 4 and mag- part of the BOSS survey (|Ross et all . l2012h . It was shown 

nitude g ~ 23, which is the objective of current and to increase by 20 to 30% the density of identified quasars, 

future cosmology projects dedicated to BAO studies in and, in particular, to effectively recover additional quasars 

the distant Universe, is a major challenge. Traditional in the redshift range 2.5 < z < 3.5. 

selections re lying on quasar colors f o r several bro a d op- Here we use this variability-based selection to measure 

tical bands (|'Schmidt fc Greed . Il983t [C room et all . 120011: the quasar luminosity function to extinction-corrected mag- 
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nitude gdorod = 22.5 and redshift z = 4 from two sets of ded- 
icated observations. For the first set of data, we took advan- 
tage of the re-observation, as part of the SDS S-III/BOSS 
survey (jEisenstein et all 1201 it iDawson et all I2012I ). of a 
14.5 deg 2 region in Stripe 82 (the SDSS Southern equa- 
torial stripe). The second s et of data used the H ectospec 
multi-object spectrograph ([Fabricant et all 120051 ) on the 
Multiple Mirror Telescope (MMT) and covered 4.7 deg 2 . 
The quasar samples being selected with only minimal color 
constrains, they are expected to be highly complete and do 
not suffer from the usual biases in their redshift distribu- 
tion induced by color selections. The selection algorithms 
of both samples are well understood, and can be applied 
to large control samples of already identified quasars in 
order to compute the completeness of the method. With 
this strategy, all corrections can be derived from the data 
themselves and do not require any modeling of quasar light 
curves or colors. 

The outline of the paper is as follows. In Sec. [5] we ex- 
plain the strategy for the selection of the targets and its 
specific application to the BOSS and the MMT observa- 
tions. In Sec. El we describe the different contributions to 
the global completeness correction for both sets of data and 
present the quasar control sample used to derive them. The 
raw and the completeness-corrected quasar number counts 
are given in Sec. HI where we also present several cross- 
checks of the results obtained. The quasar luminosity func- 
tion in g derived from these data is given in Sec. [5] 



2. Target selection 

While the basis of the target selection algorithm is the same 
for the two components of our program, it was applied with 
different thresholds to obtain the targets for the BOSS and 
the MMT observations. The BOSS component was indeed 
designed to identify a large number of quasars to a magni- 
tude limit g ~ 22.5 corresponding to the limit of BOSS 
spectroscopy at typical exposure times. The MMT data 
were primarily designed to complete the sample to fainter 
magnitudes (g ~ 23), since the telescope has a 6.5 m pri- 
mary mirror compared to the 2.5 m primary of the Sloan 
telescope. In addition, the BOSS sample was restricted to 
point sources, while the MMT sample was also used to re- 
cover quasars lying in extended sources. 



2.1. Target selection algorithm 

As a det ailed description of the variability s election can be 
found in iPalanaue-Delabrouille et al.l ([20111 ). we only sum- 
marize here the major steps of our algorithm. 

For each source, lightcurves were computed from 
the data collected by SDSS usi ng the dedi c ated S loan 
Foundation 2.5 m te lescope dGunn et all [2006). A 



mosaic CCD camera ([Grain et all 19981) ima g ed th e 
sky in the ugriz bandpasses ( Fukugita et all [l996). 
The imag ing data were proces sed through a series of 
pipelines ([Stoughton et all l2002j) which performed astro- 
metric calibration, photometric reduction and photometric 
calibration. 

The starting source list w as built from image s that re- 
sulted from the co-addition ([Annis et all I2012T ) of SDSS 



single-epoch imagesQ The completeness of the coadd cat- 
alog reaches 50% at a magnitude g — 24.6 for stars and 
g = 23.5 for galaxies. The g magnitude mentioned through- 
out the paper is computed on the coadd images, and effec- 
tively represents the mean magnitude of a variable source. 
The source morphology (point-like or extended) is also de- 
termined from these coadds. 

The lightcurves of our sources contained on average 52 
individual epochs spread over 7 years. They were used to 
compute two sets of parameters that characterize the source 
variability: 

— the x 2 s of the fit of the lightcurve by a single number 
(representing a non-varying lightcurve) in each of the 
ugriz filters (yielding five parameters, one for each fil- 
ter) 

— two para meters, an ampl i tude A and a power 7 as intro- 
duced by ISchmidt et al.l (|2010l) , that characterized the 
variability structure function V(Aiy), i.e. the change in 
magnitude Amy as a function of time lag Aty for any 
pair ij of observations, with 



V(Aty) = \Ar 



'i-j 



o? + of 



(1) 

(2) 



For each source, a neural network then combined the five 
% 2 , the power 7 (common to all filters) and the amplitudes 
A g , A r and Ai for the three filters least affected by noise and 
observational limitations (gri), to produce an estimate of 
quasar-like variability. An output j/nn of the neural network 
near designated non- varying objects, as is the case for the 
vast majority of stars, while an output near 1 indicated a 
quasar (cf . Fig. [T|) . 

This technique has been applied by BOSS for the se- 
lect ion of z > 2.2 quasars in Stripe 82, down to g ~ 

22 (|Palanque-Delabrouille et all 1201 ll ). As was clearly il- 
lustrated by this study, this approach presents the ad- 
vantage of being highly complete (the selection reached 
the unprecedented quasar completeness of 90%), even for 
a sample purity of 92%, higher than for typical selec- 
tions based on quasar colors. In addition, this variability- 
based selection was shown to overcome the drawbacks of 
color-based methods and to recover quasars near redshift 
z ~ 3 that are systematicall y missed with traditional selec- 
tions (|Richards et all l2002h . 

Here, we applied this technique to fainter magnitudes, 
with the aim of detecting quasars to g ~ 23. The magnitude 
dependence of the output of the variability neural network 
is illustrated in Fig. [TJ Requiring an output ijnn > 0.5 
selects 95% of the sample of known quasars with 18 < g < 

23 (cf. the description of this control sample in section l3~Tj) . 
Even when restricting to faint quasars with a magnitude 
g > 22, 88% of them still pass this threshold. 

In addition to the variability-based selection, a loose 
color constraint was used to reject a region of color-space 
mostly populated by stars, in order to reduce the fraction 
of stellar contaminants in the target list. The cut consisted 
in requiring C3 < a — ci/3 (cf. Fig. [5]), where the level a was 
set to 1.0 or 0.6 depending on the target morphology (see 
Sec. 12.21 and 12.31 for details), and c\ and C3 are defined as 



We used t he Catalog Archive Server (CAS) interface 
(http://casjobs.sdss.org) to recover both the Stripe 82 coadd 
and the single epoch information. 
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Fig. 1: Output of the variability Neural Network as a func- 
tion of g magnitude for a sample of known stars (small black 
dots near yNN = 0) and for known quasars (larger red dots 
at y NN ~ 1). 



mUml rfl999h by: 

ci = Q.95(u-g) +0.31(5 - r) + 0.11(r-«), 

c 3 = -0.39(u - 5 ) + 0.79( ff - r) + 0.47(r - i) . (3) 

About 100% of z < 2.2 and 98% of z > 2.2 known quasars 
(resp. 95% and 93%) pass the condition with a — 1.0 (resp. 
a = 0.6). 




-0.5 0.5 1 1.5 2 2.5 3 3.5 4 



Fig. 2: Locus of stars (upper thin blue contours), z < 2.2 
quasars (lower left thick green contours) and z > 2.2 
quasars (lower right thick red contours) in the C3 vs. c\ 
color-color plane. The upper solid line corresponds to the 
color cut C3 < 1.0 — Ci/3 and the lower dashed line to 
c 3 < 0.6 -a/3 (cf. Eq.E}. 



2.2. The BOSS sample 

The BOSS fields dedicated to this high density quasar sur- 
vey were part of Stripe 82 and located at 317° < aj2ooo < 
330° and 0° < S J2 ooo < 1-25°. 



The starting source list consisted of all point sources 
in this area passing the usual BOSS quality cuts (on pho- 
tometry, pixel satur ation, sourc e blend ing etc.) as described 
in the Appendix of iRoss et al.l (|2012l ). as well as the loose 
color cut (with a = 1.0) mentioned in the previous section, 
yielding about 5200 objects per deg 2 . 

All objects passing yNN > 0.5 were selected as targets. 
The resulting list contained ~ 270 targets per deg 2 . We 
removed from this list all tar gets t hat had already been 
observed by BOSS (|Ahn et al.U2012D . This reduced the list 
to ~ 240 objects per deg 2 . 

In case of fiber collisions during the tiling procedure, 
priorities were set on the targets according to their mag- 
nitude, giving higher priority to the brightest objects since 
these are more likely to obtain an accurate identification. 
An identical maximal priority was attributed to all targets 
with g < 22.7 (not corrected for extinction as the quantity 
of interest here is the actual observed flux of the object), 
and lesser priorities for fainter magnitudes. 

Seven partially overlapping BOSS half-plates were allo- 
cated to the project (numbered 5141 to 5147 from low to 
high aj2ooo)> covering a total of 14.5 deg 2 . The other plate 
halves were used for the standard BOSS survey and are not 
included in this analysis. Plate 5141 was exposed for over 
2 hours (ie. twice the normal time), without significant in- 
crease in the number of quasars identified. The subsequent 
plates were thus exposed for 1 hour. All seven plates were 
observed during July and September 2011. 

2.3. The MMT sample 

For the sake of completeness, the starting list for the MMT 
sample was built from all SDSS sources in the area of the 
previous BOSS sample, whether point-like or extended and 
whether or not they passed the standard BOSS quality cuts, 
and passing the loose color cut with a = 1.0 described in 
Sec. 12. 11 This initial list consisted of ~ 11300 deg -2 objects, 
60% of which were resolved (ie. "extended", like galaxies), 
and 40% of which were unresolved (ie. "point-like", like 
stars). As for the BOSS sample, all objects that have al- 
ready obtained spectra with SDSS-III/BOSS were removed 
from the target list. In addition, we rejected targets with 
g < 22 that were simultaneously selected for the BOSS 
sample, since the redshift determination efficiency of BOSS 
is close to unity at least to that magnitude. 

Most quasars are expected to be significantly more lu- 
minous than their host galaxy, making them appear point- 
like. We thus expect less quasars in the "extended" sub- 
sample. The fraction of random objects artificially drifting 
into the region of large yNN values being independent of 
the source morphology, however, the ratio of the number of 
quasars to the total number of selected targets is therefore 
smaller for extended sources. We thus applied stricter cuts 
to the extended sub-sample, which is affected by a higher 
contamination, and we gave a higher priority to the MMT 
point-like targets. 

The highest priority (priority A) was given to point-like 
targets with observed magnitude g < 23.0, yNN > 0.5 and 
C3 < 1.0 — 0.33ci, leading to ~ 130 deg -2 targets. 

Priority B was given to extended objects with g < 23.0, 
y N N > 0.8 and c 3 < 0.6 - 0.33 x ci. This led to 140 deg" 2 
additional targets. 
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Priority C (respectively D) was given to objects with 
23.0 < g < 23.2 passing the same conditions as the priority 
A (resp. B) targets, leading to 40 deg -2 (resp. 30 deg~ 2 ) 
additional targets. 

A total of 6 non-overlapping circular fields of 1 deg in 
diameter were observed with MMT/Hectospec in the last 
trimester of 2011, with exposure times of 150 minutes each. 
The fields were located within the coverage of the BOSS 
sample of Sec. 12.21 in the area farthest from the Galactic 
plane, which is least contaminated by Galactic stars. The 
region 326.5° < cfj2ooo < 328° was avoided because of a 
higher extinction on average (A g ~ 0.4) than for the other 
MMT fields (A g ~ 0.2 for 322° < a J2 ooo < 326.5° and 
A g ~ 0.3 for 328° < a. rannn < 330°), where t he extinctions 
come from the maps of lSchlegel etafl ffl998h . 

Because of weather conditions, only five MMT fields (all 
but the field with lowest Right Ascencion) produced us- 
able data. This resulted in a total sky coverage of 3.92 deg 2 . 

Figure E] illustrates the locations of the BOSS and MMT 
fields dedicated to the present study. 
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Fig. 4: Fraction of known quasars identified as point-like in 
co-added images, as a function of magnitude, for several 
bins in redshift (rising redshift from bottom to top). The 
asymptotic red histogram is the fraction of stars that are 
identified as point-like. 




326 328 330 

R.A. (J2000) degrees 

Fig. 3: Footprint of the BOSS half-plates (black contours) 
and of the MMT fields (filled red disks) dedicated to the 
present study. BOSS plates were attributed the numbers 
5141 to 5147 from left to right. MMT fields were labeled 
to 5 with increasing Right Ascension. 



frame, pointing towards a physical effect that cannot be ex- 
plained with the previous arguments. Even the curves for 
2.0 < z < 3.0 quasars are significantly lower than that ob- 
tained for stars, our control sample of point-like objects. 
There is a clear step at a redshift z ~ 0.8, where 60 to 80% 
of the quasars, even bright (g < 20) ones, appear extended. 

These results are consistent with the host galaxy becom- 
ing detectable in the co-added images, making the quasar 
appear extended. Fig. @] thus indicates a statistical detec- 
tion of the host galaxy of quasars, even at redshifts as large 
as z ~ 2. The brightest (g < 19) quasars, however, still 
sufficiently outshine their host galaxy to remain point-like, 
except in the lowest redshift bin where the host galaxy is 
resolved. 



3. Completeness corrections 



2.4. Impact of source morphology 

While most quasars, whatever their redshift, appear point- 
like on single-epoch SDSS images, this is no longer the case 
on co-added frames, which result from the superposition, at 
the image level, of about 56 Stripe 82 scans, thus reaching 
a depth about 2 magnitudes fainter than individual scans. 
In the BOSS sample, we thus rejected a large fraction of 
low-redshift quasars by only considering point-like objects. 
Fig. [4] shows the fraction of quasars classified as point-like 
on the co-added images, for several bins in redshift. 

Fig. @] also displays the fraction of point-like objects 
on a single-epoch image of excellent seeing (most of these 
objects are indeed stars) that are classified as stars on the 
co-added frame. While all stars with magnitudes g < 22 ap- 
pear point-like in the deep frame, up to 16% of them appear 
extended at g = 23, possibly due to small misalignments in 
the co-addition procedure that affect the object morphology 
at faint magnitudes. High redshift (z > 3) quasars follow 
the same trend as stars, therefore suffering from the same 
technical drawbacks. In contrast, as the redshift decreases, 
more and more quasars appear extended on the co-added 



To estimate quasar densities and the luminosity function 
from raw quasar counts, the data have to be corrected for 
biases introduced by the cuts applied to select the targets 
and for instrument-related losses. The analysis-related com- 
pleteness corrections are computed from control samples of 
known quasars, while the instrument-related ones are com- 
puted from the data. 

The different contributions to the completeness correc- 
tions are detailed in the following sections, and a summary 
of the magnitude dependence of the corrections for BOSS 
and the MMT is illustrated in Fig. [5] 

3.1. Control sample of known quasars 

The analysis-related completeness corrections are deter- 
mined using a list of 19215 spectroscopically confirmed 
quasars in Stripe 82 obtai ned from the 2dF quasar catalog 
(|2QZ: Croom et alll200i. the 2dF-SDSS LRG and Quasar 
Survey 2SLAQ (ICroom et all 120091) the SDSS-DR7 spec- 
troscopic datab ase (lAbazaiian et alll2009r) . the SDSS-DR7 
quasar catalog (I Schneider et all 1201 Oft and BOSS observa- 
tions up to July 2011 (|Ahn et all 120121: iParis et alll2012h . 
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Fig. 5: Magnitude dependence of the contributions to the total completeness correction for the different samples. Upper 
left: the new quasars identified in the BOSS sample. Upper right: the known quasars that were intentionally removed 
from the target lists. Lower left (resp. right): the quasars identified in the MMT sample and that appeared point-like 
(resp. extended) in the coadded SDSS frame. These corrections are computed from all quasars contained in either sample. 



As shown in Fig[6j these quasars have redshifts in the range 
< z < 5 and g magnitudes in the range 17 < g < 23 
(Galactic extinction-corrected) with 212 quasars having a 
magnitude fainter than g = 22.5. The irregular shape of the 
distributions results from the use of several surveys with 
different redshift goals and selection algorithms. At z > 2, 
most of the quasars c ome from color-selection with the stan- 



dard BOSS survey (lYeche etHl 120101: iKirkpatrick et all . 
120111: IBovv et all . 120111: iRoss et all . l2012h . 

Approximately 2% of the objects could not be matched 
to any so-called primary target in the SDSS database, re- 
ducing the sample to 18910. Lightcurves and variability pa- 
rameters were computed for these quasars using the proce- 
dure described in Sec. [2] 



3.2. Analysis-related completeness corrections 

The analysis-related corrections arise from two contribu- 
tions that affect the BOSS and the MMT samples differ- 
ently: a quality factor e qua i and a selection factor e so i. These 
corrections depend on the source morphology (point-like or 
extended), the redshift z and the magnitude g uncorrected 
for extinction (since we are here sensitive to the observed 



and not to the intrinsic quasar magnitude). For the sake 
of clarity, the magnitude corrected for Galactic extinction 
will henceforth be denoted f/dered- 



Quality completeness e qua i(g): 

This contribution i s related to the q uality cuts described 
in the Appendix of IRoss et al.l (|2012f ) that were applied to 
obtain the initial list of the sources from which the selection 
was made. It only affects BOSS data as no quality cut was 
applied to build the list of sources for the MMT sample 
(see sections 12.21 and |2~3|) . The mean quality factor over the 
BOSS data of the present study is 0.89, with a plateau at 
0.90 for all bright sources with g < 22, dropping to 0.70 at 
g = 23 (see blue curve in upper two plots of figure [5]) • 



Target selection completeness 

esei(<7, z , source morphology) : 

This contribution is related to the variability-based and 
color-based selection algorithms, as described in Sec. [2] It 
depends on the telescope that the target list is designed for 
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Fig. 6: Distributions in redshift (top plot) and extinction- 
corrected magnitude (bottom plot) for the sample of 
quasars from previous quasar surveys covering Stripe 82. 



only because of the different thresholds that were set on the 
selection variables. 

BOSS targets were selected from an initial list limited 
to point-sources. To take into account the completeness loss 
due to the non-consideration of extended sources, e sc iBOSS 
was computed from the ratio, in the control sample, of the 
point-like quasars that passed the selection cuts to the total 
number of quasars. This yielded the correction table shown 
in Fig. [7] as a function of magnitude g and redshift z. The 
mean completeness correction for the selection step, over 
the BOSS sample, is 0.78. 

For the MMT data, the completeness correction for 
the selection step was computed separately for point-like 
sources and for extended ones. For each source morphol- 
ogy, £ so immt(3, z) was determined from the control sample 
as the fraction of quasars of a given morphology that passed 
the relevant selection criteria (cf. Sec. I2.3[) . The mean se- 
lection completeness correction over the MMT sample of 
quasars is 0.76 (0.82 for the point-like sources and 0.57 for 
the extended ones). 

3.3. Instrument-related completeness corrections 

These corrections come from two contributions: a tiling fac- 
tor etiiing quantifying whether a target could indeed have 
been observed, and a spectrograph factor e spe ctro quanti- 



Function using dedicated data from BOSS and MMT 




Observed g magnitude 



Fig. 7: Selection completeness e se iBOSs(ffi z ) for the quasars 
in the BOSS program. The significant drop at low redshift 
is due to the fact that many quasars with z < 1 appear 
as extended sources in coadded frames and were thus not 
considered in the initial source list. We also observe the 
expected efficiency drop at faint magnitudes. Bins are left 
unfilled (white) are bereft of quasars. 



fying the probability of obtaining a secure identification 
and redshift for a given target. These depend on the source 
morphology (point-like or extended) through the target pri- 
ority, and on the observed magnitude g. 

Tiling completeness e t m n g (priority): 

This contribution is related to the tiling proce- 
dure (jDawson et all l2012t) . taking into account the den- 
sity of allocated fibers and the target priorities in case of 
fiber collisions. This correction just represents the fraction 
of targets that were assigned fibers. It entirely depends on 
the instrumental configuration and thus must be computed 
specifically for BOSS and the MMT. As explained in Seed 
the target priority varies with magnitude for BOSS targets, 
and with both magnitude and source morphology (point- 
like or extended) for MMT targets. In the case of BOSS, 
however, the fraction of spectra that can be correctly iden- 
tified drops significantly beyond g ~ 22.5 and only the first 
priority bin was used for the present analysis (cf. Sec. I4.3[) . 
The mean tiling completeness correction is 0.95 for the 
BOSS data and 0.84 for the MMT data. Corrections per 
priority bin are given in table [T] 



Project 


BOSS 


MMT 


MMT 


MMT 


MMT 


(Priority) 


A 


B 


C 


D 


^tiling 


0.95 


0.90 


0.74 


0.68 


0.53 



Table 1: Tiling completeness correction for BOSS (unique 
priority level) and the MMT (priorities A through D). 



Spectrograph completeness e sp cctro(5): 

Some spectra did not produce a reliable identification of the 
source, either because the extraction procedure had failed 
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(hereafter "Bad" spectra) or because the spectrum had too 
low a signal-to-noise ratio for adequate identification (here- 
after "?" ) . Fig. [5] illustrates the rate of the "bad" spectra 
as a function of magnitude g for the BOSS (filled black 
dots) and the MMT (open red circles) spectra. The over- 
laid dashed curves are fits to the data. The fraction of "bad" 
spectra reaches 15% at g = 23.15 for the MMT and g = 22.9 
for BOSS. The upper blue dot-dashed curve illustrates the 
total loss 1 — e S pcctro(9), similar for both instruments, when 
considering both the "bad" and the "?" spectra. This total 
fraction of inconclusive spectra is negligible for g < 22, it 
reaches 15% at g = 22.8, and 30% at g = 23.1. 




Fig. 8: Fraction of inconclusive spectra (declared so after 
visual inspection) as a function of observed g magnitude. 
The curves are fits to the data. In our analysis, we use 
the upper blue curve to correct for the total fraction of 
inconclusive spectra. 



4. Results 

BOSS data were taken on th e dedic ated Sloan Foundation 
2.5 m tele scope (iGunn et all 120061 ) using the BOSS spec- 
trograph (|Smee et all l2012f) . The BOSS spe ctra were re 
duced with the stan d ard B OSS pipe line dBolton et al 



2012 1 ) : iDawson et all (|2012h . see also Bolton & Schlegcl 



20101) ) which also provides an automated determination 



of the target classification and redshift. A manual inspec- 
tion was performed on all the spectra of this program, as 
for al l quasar targets of the main BOSS survey ()Paris et all 
12012ft . The MMT spect ra were reduced with a customized 
pipeline based heavily on the E-SPECROAD packageQ All 
the spectra were checked visually to produce final identifi- 
cations and redshifts. The spectra were classified as "QSO" 
for secure quasars with reliable redshift, "QSO?" for se- 
cure quasars but uncertain redshift, "Star", "Galaxy" or 
"Inconclusive" . The latter case encompasses the "bad" and 
the "?" spectra of the previous section. In our analysis, we 
use all spectra that were identified as "QSO" or "QSO?" . 



4.1. Raw counts for the BOSS and MMT samples 

We identified 1179 new quasars with BOSS (hereafter re- 
ferred to as "New BOSS" or simply "BOSS") and 262 ones 
with the MMT. To these, we add the 436 previously iden- 
tified quasars in our area of interest that were explicitly 
removed from the target lists and hereafter referred to as 
"Known" . 

Fig. IH1 shows the distributions of the redshift and ob- 
served magnitude in the g band for our sample of confirmed 
quasars. While the redshift distributions of the BOSS and 
MMT samples are similar in shape, the MMT sample 
reaches more than half a magnitude deeper than the BOSS 
sample. The average magnitude of each of the two samples 
is (j)mmt = 22.2 and (s}boss = 21.5. The sample of pre- 
viously known quasars that were not included in the target 
lists have a mean magnitude (g) Known = 21.0. 




22 23 24 

Observed g magnitude 




http:/ /iparrizar. mnstate.edu/ juan/research/ESPECROAD/ 



Fig. 9: Quasar number counts per deg 2 as a function of ob- 
served g magnitude (top) and redshift (bottom) . The black 
(resp. red) line is for quasars measured by BOSS (resp. 
MMT). The BOSS histograms include the already identi- 
fied quasars in the area that were removed from the target 
list. 



As explained in Sec. H21 the MMT targets at g < 22 
only consisted in the selected objects that were not included 
in the BOSS sample, i.e., either that did not pass the quality 
criteria or that appeared as extended in the source catalog. 
This led to 69 quasars with g < 22 in the MMT sample. For 
g > 22, the MMT sample included 116 quasars in common 
with BOSS, and 77 additional ones. 
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Table [5] summarizes the raw quasar counts in the BOSS 
and MMT samples. 



Sample 


< 22 


Observed g 
22 - 22.5 


magnitude 
22.5 - 23 


> 23 


Total 


BOSS 


692 


231 


165 


91 


1179 


MMT 


69 


88 


73 


32 


262 


Known 


381 


51 


4 





436 



Table 2: Raw number counts for the different samples in 
several g magnitude bins. 



4.2. Identification cross-checks 

We have intentionally observed some targets with both in- 
struments in the g > 22 magnitude regime where we ex- 
pect the BOSS identification to become less secure. Table[3] 
summarizes the cross-identification of the 344 common tar- 
gets, among which 116 quasars identified as such from both 
instruments. There are very few changes in identification. 
The most notable feature is that, as expected, more faint 
quasars could be identified and given a secure redshift from 
MMT than from BOSS: 19 "QSO?" BOSS targets were con- 
firmed as "QSO" with MMT, as well as 9 targets classified 
as "Bad" in BOSS. It is noteworthy that even at these faint 
magnitudes, there were almost no false quasar detections 
in BOSS: only 1 "QSO?" BOSS target was identified as a 
star from the MMT spectrum; all others were confirmed as 
quasars. 



MMT\BOSS 


QSO 


QSO? 


Star 


Galaxy 


Unknown 


QSO 


94 


19 


3 




9 


QSO? 


2 


1 


3 




2 


Star 




1 


115 


1 


37 


Galaxy 






6 


6 




Unknown 


2 


3 


19 




21 



Table 3: Cross-identification between BOSS and MMT 
spectra for the 344 targets with g > 22 that were observed 
by both telescopes. 



At brighter magnitudes (22 < g < 22.5), table U sug- 
gests excellent consistency in the identification by either 
telescope. Out of a total of 66 targets for which a secure 
identification is available from the MMT spectra, 59 were 
correctly identified by BOSS, 5 had too low S/N, and 2 
were misidcntified. There are no false quasar detection. 



MMT\BOSS 


QSO 


QSO? 


Star 


Galaxy 


Unknown 


QSO 


38 


5 


2 




4 


QSO? 










1 


Star 






15 






Galaxy 








1 




Unknown 






2 







Table 4: Cross-identification between BOSS and MMT 
spectra for the 68 targets with 22 < g < 22.5 that were 
observed by both telescopes. 



From the set of targets with an identification in both 
samples (i.e., not considering those declared "unknown" in 
either case) , and assuming the true identification to be the 
one from the MMT spectrum, we can estimate the iden- 
tification reliability with BOSS to be of order 237/251 = 
94 ± 2% at g > 22 and of order 59/61=97 ± 2% over the 
magnitude range 22 < g < 22.5. 

4.3. Corrected counts 

To study quasar number counts, we use the quasars iden- 
tified in the BOSS sample to g < 22.5 and in the MMT 
sample to g < 23.0. The depth of the observations is in- 
sufficient to use identified quasars with g > 23.0. All fields 
suffer from relatively high Galactic extinction: (A g ) = 0.27 
(rms of 0.10) ranging from 0.15 to 0.60 in the BOSS fields 
and (A g ) = 0.25 (rms of 0.04), ranging from 0.15 to 0.35, 
in the MMT fields. The previous limitations on observed 
magnitude therefore result in the following effective limits 
on the extinction-corrected magnitude: <?dercd < 22.35 for 
BOSS and g deIed < 22.70 for the MMT. 

To compute completeness-corrected quasar counts, we 
defined three exclusive areas, based on observed g magni- 
tude and sky coverage: 

- g < 22.0 : use of the BOSS sample over the entire 
14.5 deg 2 area (691 BOSS and 346 known, for a total of 
1037 quasars) 

- 22 < g < 22.5 : use of the MMT sample over its 3.9 deg 2 
coverage (88 quasars) and of the BOSS sample for the 
remaining 10.6 deg 2 area (169 quasars), completed by 
47 known quasars 

- 22.5 < g < 23.0 use of the MMT sample over its 
3.9 deg 2 coverage (105 MMT and 3 known quasars) 

This division ensures, for each magnitude range, the use 
of the sample with best redshift reliability (cf. Sec. 14. 2\i 
and maximal statistical significance. The data that did not 
explicitly enter the computation of the corrected counts or 
quasar luminosity function were used for cross-checks (cf. 
Sees. EM and EM- 

The data from the BOSS sample were corrected for se- 
lection (e se i) and quality (e qua i) incompleteness. In addition, 
tiling (etiimg) and spectrograph (e spe ctro) completeness cor- 
rections were applied to the quasars that were identified 
from this deep program. The previously identified quasars 
in the area were only corrected by e se i and e qua i since they 
were removed from the list prior to the tiling procedure. 

The data from the MMT sample were corrected for 
spectrograph (e spcc tro) incompleteness, and for the relevant 
etmng and e SG i that, for the MMT sample, depended on the 
target morphology. 

The completeness-corrected number of quasars is thus 
computed from the following equation: 

x 1 

Nqso = 2^ 

Ar ^scl ^qual ^tiling ^spectro 

JVBOSS 

+ £ — — 

AT ^scl ^qual 

Known 

+ £ 



JV. 



point — like sel 



point— like 



^tiling ^spectro 
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E 



i 



^extended 
fc sel 



^tiling ^spectro 



(4) 



where the completeness corrections are 2D functions of the 
quasar magnitude g and redshift z. 

Figure [TO] illustrates the completeness-corrected quasar 
number counts as a function of redshift f or 3 magnitude 
limits : gdcrod < 22 as in the BOSS survey (|Schlegel et al l. 
2009), 5dcrcd < 22.5 as required for the eBOSS projedQ, 
and <?dered < 23 as required for BigBOSS__. As explained 
above, however, the last set of histograms are only com- 
plete to .gdered < 22.70 and thus represent lower limits for 
ffdered < 23. On average over the quasars of this project, 
total completeness corrections are at the level of 80% to 
g < 20, and drop smoothly to 50% at g ~ 22.5. Table[S]sum- 
marizes the total quasar number counts for various ranges 
in corrected g magnitude and redshift. 



Redshift 

z 



< 1 
[1 - 2.15] 
> 2.15 
all z 



Extinction-corrected magnitude (jdered 
< 22 < 22.25 < 22.5 < 22.75 < 23 



27(3) 
75(3) 
34(2) 
136(5 ) 



29(3) 
87(3) 
37(2) 
153(5) 



30(4) 
99(4) 
47(3) 
175(6) 



30(4) 
113(5) 
49(3) 
191(7) 



32(4) 
119(6) 
53(3) 
205(8) 



Table 5: Number counts (quasars per deg 2 ) for several 
ranges in extinction-corrected magnitude gdered and red- 
shift z. The statistical uncertainty on the last significant 
digit is indicated in parentheses. Number counts are com- 
plete to <7dcrcd < 22.5, and only indicate lower limits at 
fainter magnitudes. 



4.4. Counting cross-checks 

Because we had overlapping data in terms of magnitude 
range or sky coverage from two different programs, sev- 
eral counting cross-checks can be performed. Every single 
quasar observed was therefore used to compute the lumi- 
nosity function or for cross-checks, and some for both. 



Density comparison in 22 < g < 22.5 

The magnitude range 22 < g < 22.5 was accessible to both 
the MMT and the BOSS samples. Because of the differ- 
ent constraints on source morphology and quality, however, 
and since the MMT zone only covers a portion of the BOSS 
zone, the quasars observed in the two cases were not the 
same. The completeness-corrected quasar number counts 
derived from either sample over 22 < g < 22.5 are given in 
table El The last column indicates the number of identified 
quasars used to make the measurement. The quasar densi- 
ties derived from BOSS or the MMT over this magnitude 
range are in excellent agreement. 



Survey 


QSO density 
(deg" 2 ) 


Raw number 
of quasars 


BOSS over MMT zone 


39.0 ± 5 


72 


BOSS complete sample 


40.2 ± 3 


276 


MMT 


39.1 ± 4 


100 



Table 6: Efficiency-corrected number counts for the BOSS 
and MMT samples (including, in both cases, the BOSS 
standard QSOs that were selected but not re-observed) over 
the magnitude range 22 < g < 22.5 accessible to both sam- 
ples. 




01234S1234912345 

Redshift 



Fig. 10: Quasar number counts per deg 2 as a function of 
redshift (Az = 0.2). Blue is for the new sample of quasars 
identified with our BOSS deep program, green for the pre- 
viously known BOSS sample, red is for the MMT sample 
and black for the total. The blue, green and red curves use 
mutually exclusive samples and correspond to the zones 
defined in Sec. 14.31 



3 http:/ /lamwws. oamp.fr/cosmowiki/Project_eBoss 

4 http://bigboss.lbl.gov 



Cross-check of e qua i using MMT point-like targets at g < 22 

The sample of quasars selected for the MMT observations 
from point-like sources in the magnitude range g < 22 
should allow recovery of the quasars that were not se- 
lected in the BOSS target list because of the constraints on 
source quality. They can therefore be used to verify the es- 
timated e qua i. This concerned N = 22 MMT quasars, with 
a mean tiling correction (corresponding to priority level A) 
^tiling = 0.903. Given a survey area S = 3.9 deg 2 , the den- 
sity of point-like sources observed with MMT at g < 22 is 
therefore N/e mm /S = 6.2 ± 1.3 deg~ 2 . 

In the same magnitude range, the BOSS sample con- 
sisted of N = 691 quasars, with a mean tiling correc- 
tion eyiing = 0.945 and a mean source quality complete- 
ness correction e qua i = 0.897. Given a survey area S = 
14.5 deg 2 , the estimated density of quasars not included 
in the sample because of the quality constraint is thus 
N/etiUng/equal x ( l ~ e quai) / 'S = 5.8 ± 0.2 deg -2 , in agree- 
ment within la with what was estimated from the MMT 
sample. 

Cross-check of e sc i using MMT extended targets at g < 22.5 

Finally, the quasars selected for the MMT observations 
from extended sources at g < 22.5 should allow recovery 
of the quasars that were not selected in BOSS because of 
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constraints on the source morphology. This can be checked 
by comparing the density of extended quasars in the MMT 
survey to the number of quasars that were not selected in 
BOSS for this same reason, estimated from the observed 
number of point-like quasars and the morphology-part of 
the selection completeness correction. 



o 

a 10 



6 — 
4 — 

2 — 



-j Extended targets in MMT 

Morphology correction (BOSS) 



0.5 



1.5 



2.5 



3.5 
Z 



Fig. 11: Comparison of the number of "extended" quasars 
estimated from BOSS data and completeness corrections 
(blue curve) or observed in the MMT sample (purple 
points). 



Quantitatively, we compute the following two quanti- 
ties, in the g < 22.5 magnitude range. On the one hand, 
the density of extended quasars identified in MMT: 



1 



E 



5mmt e' , 

/vrcxtcndcd SCI 
JV MMT 



extended 



^tiling ^spectro 



and on the other hand, the density of nnselected extended 
qnasars in BOSS: 



1 



£boss J-^ ^tiling ^spectro ^qual 
-"BOSS 



1 



Si 



BOSS 



iV Kn 



^qual 



"morphology 



^morphology 
\ 



l 



where Smmt an d <5boss ar e respectively the areas of the 
MMT and the BOSS programs, and e m orphoiogy is the frac- 
tion of the targets that are point-like (as a function of mag- 
nitude and redshift). This correction is the part of e se i that 
does not include the effect of the target selection based on 
its color and time variability. These two densities are illus- 
trated in Fig. 111! They are clearly in agreement. 



5. Luminosity function in g 

We compute the quasar luminosity function (LF) from the 
corrected number counts derived above, and considering 
our completeness limit at (?rjered < 22.5. The distance mod- 
ulus cLm{z) is computed using the standar d flat ACDM 
model with the cosmological parameters of lLarson et al.l 
(|2011h : n M = 0.267, n A = 0.734 and h = 0.71. 



5.1. K-corrections 

Selection for this survey was performed in the <?-band, and 
for the majority of the data this band provides the high- 
est S/N. We define the if-cor rection in terms of the ob- 
served g magnitude and follow [Croom et aD (|2009f ) (here- 
after C09) in applying the correction relative to z = 2, 
whic h is near the median r edshift of our quasar sample (see 
also iRichards et al.l (|2009h ). The absolute magnitude nor- 
malized to z = 2 is given by: 



M g {z = 2) = g c - d M {z) - \K{z) - K{z = 2)] 



(5) 



Hereafter we will use M g as a shorthand for the redshift- 
corrected M g (z = 2). The if-correction as a function of 
redshift is de rived from model qu asar spectra in a simi- 
lar fashion to Ric hards et al.l (|2009f) . The quasar model in- 
cludes a broken power law continuum with ay = —0.5 at 
A > HOOA and a y = -1.5 at A < HOOl (|Telfer et all 
2002). Strong quasar emission lines are included, where 
the equivalent width is a function o f luminos i ty ac cording 
to the well-known Baldwin Effect (Bal dwin! . Il977l ): thus, 
the X-correction is a function not only of redshift but 
also luminosity (or equivalently, observed g-magnitude) . 
The model also includes Fe em ission using the template of 
IVestergaard and Wilked (120011) and Lyman-q forest absorp - 
tion using the prescription o f lWorseck fe Prochaskal (|2010h . 
The forest model is particularly relevant here, as for z ^ 2.5 
the g-band if -correction necessarily includes a component 
due to forest absorption; our if -correction accounts for the 
mean value but for high redshift objects the uncertainty in 
individual if-corrections is increased by line-of-sight fluc- 
tuations in the amount of forest absorption within the g- 
band. The models used to derive the if -corrections will be 
described in fuller detail in McGreer et al. 2012 (in prep). 
In general, as shown in Fig. 1121 the values are very sim- 
ilar to those used by C09 (e.g., their Fig. 1), but are ex- 
tended to z — 4. At z ~ 2 — 3, the luminosity-dependent 
if-correction we use introduces a differential value of ~ 0.25 
mag across the luminosity range of quasars in our sample, 
as the Lyman-a and C IV lines are within the g-band and 
contribute substantially to the flux within the bandpass. 



5.2. Luminosity function model 

We define eight redshift bins: the first five are 
the same redshift intervals as in C09, with limits 
0.68,1.06,1.44,1.82,2.2,2.6; the last three are specific to 
our analysis and have the limits 2.6,3.0,3.5,4.0. We calcu- 
late the bin ned LF using the m odel-weighted estimator $ 
suggested bv lMivaii et all (|200lh . which presents the advan- 
tage of not having to assume a uniform distribution across 
each bin, unlike 1/V estimators. Instead, it models the un- 
binned LF data and uses it to correct for the variation of 
the LF and for the completeness within each bin, which is 
here particularly critical at the faint end of the LF where 
the latter is incompletely sampled. This estimator gives the 
binned LF as 



odd 



{M 9i ,Zi) 



N, 



obs 



odel 



(6) 



where M gi and Zi are, respectively, the absolute magni- 
tude and the redshift at the center of bin i, § model is the 
model LF estimated at the center of the bin, ]\l] nodel is 
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B 1-6 

0) 



K(z) for this work 

K(z) = -2.5(1+vJlog (1«) 



/ 



- 



LF measurements are all systematically below the best-fit 
C09 model and this trend is confirmed by our data in the 
same redshift bin, confirming the hypothesis that the fea- 
ture is real and not related to a redshift limitation of the 
sample since ours extends significantly beyond z — 2.6. 
This disagreement was already noted by C09 who observed 
a significant improvement of the fit when reducing the up- 
per limit in redshift to z < 2.1. 

5.3, Luminosity function fit 

The data from our analysis and from C09 are indepen- 
dent. We therefore combine them to constrain the luminos- 
ity function in larger redshift and magnitude ranges than 
either data set alone. 

Fitting the model defined in Eqs. [7] and \E\ to the com- 
bined LF does not significantly change the result from the 
best-fit model of C09. In particular, Eqs. [7] and [5] do not 
provide sufficient freedom to solve the discrepancies men- 
tioned above. We therefore defined a new fitting model, 
based on the same equations as before but allowing the 
redshift evolution parameters (ki and k 2 ) and the model 
slopes (a and (3) to be different on either side of a pivot 
redshift z p — 2.2. The model is thus described by Eq. [7J 
the number of quasars estimated from the model in the where a and (3 now have subscripts I for z < z p and h for 



-0.5 



1.5 



2.5 



3.5 



Fig. 12: K-corrections as a function of redshift. The spread 
illustrates the variation over the luminosity range of the 
quasars in our sample. The red curve shows for comparison 
the parameterization K (z) = —2.5(1 + v a ) log 10 (l + z) in- 
troduced in C09, with v a = —0.5 and normalized to z = 2. 



bin (accounting for completeness corrections) and N° is 
the observed number of quasars in the bin. Systematics on 
the LF measurements are computed by estimating the un- 
certainty on the completeness corrections of Sec. [3] using 
Monte Carlo simulations in which the selected and total 
number of quasars follow Poisson distributions centered on 
the observed means. 

We assume the LF to be appropria tely described by a 
standard double power-law of the form (jBovle et all [2000) : 



z > z p , and an evolution in M* characterized by 
M*(z) = M* g {z p ) - 2.5[fci(* - z p ) + k 2 (z- z p ) 2 } , 



(9) 



10 0.4(a+l)(M 9 -M*) + 10 0.4(/3+l)(M g -M«) 



(7) 



where <& is the quasar comoving space density. Following the 
pure luminosity evolution model of C09, a redshift depen- 
dence is introduced through an evolution in M* described 
by 



M*Jz) = M.*(0) - 2.5(fci2 + k 2 z 2 



(8) 



In a first step, we impose, for our model $ mode ^ the pa- 
rameters fitted by C09 for the 0.4 < z < 2.6 redshift range 
and given in their table 2, ie. M*(0) = -22.17, ki = 1.46 
and k 2 = —0.328. Figure Q1J] illustrates our LF measure- 
ments (black circles) where the error bars include both our 
statistical and systematic uncertainties added in quadra- 
ture. The blue open circles are the LF measurements from 
C09 and the blue dashed curve is their best-fit $ mode/ . 

The data presented in this work from dedicated BOSS 
and MMT surveys cover the redshift range from z — 0.68 
to z = 4, allowing 3 additional redshift bins beyond those 
of C09. Over the 0.68 < z < 2.6 redshift range common to 
both studies, our data and the C09 points are in good agree- 
ment. The extrapolation of the C09 model to z = 4 is also 
in reasonable agreement with our data points. Two discrep- 
ancies, however, are worth being noted. In the 3.5 < z < 4.0 
redshift bin, our measured LF points are all above the ex- 
trapolated curve; the uncertainties on the measured LF, 
however, are very large and the trend needs to be confirmed 
with improved statistics. In the 2.2 < z < 2.6 bin, the C09 



where k\ and k 2 are again to be considered separately for 
low (subscript I) and high (subscript h) redshifts w.r.t. z p . 

The best-fit parameters are reported in Table [Jj and the 
resulting luminosity function is illustrated in Fig.[T2]as the 
red curve. The best fit has a \ 2 of 324 for 162 degrees of 
freedom. For the sake of comparison, we also provide the 
parameters of C09, translated into the parameters defined 
in Eq. We provide, in Appendix A, the measured lumi- 
nosity function values and associated uncertainties, where 
the measurements were corrected using the estimator de- 
scribed in Sec. 15.21 and considering the best-fit model. 

In the redshift bins from z = 0.68 to z < 2.6 where the 
data from C09 are the most constraining, the best-fit model 
shows negligible difference with that of C09. It starts to de- 
viate at redshifts z > 3.0, with a significantly flatter faint 
luminosity slope (smaller \/3\) than observed at low redshift. 
Further disc ussions on t he quas ar luminosity function can 
be found in iRoss et ail (|2012bf) . where the BOSS data of 
this analysis are combined with the full data set from the 
DR9 release of the SDSS-III/BOSS survey. The quasar se- 
lection there requires models of quasar colors to determine 
the selection completeness, but the statistics are increased 
by over an order of magnitude, thus allowing better con- 
strains on the LF model in the high redshift bins. 

Using our best-fit parameters of the quasar luminosity 
function, we provide, in Tabled an estimate of the number 
of quasars in the redshift-magnitude plane for an hypo- 
thetical survey covering 10,000 deg 2 . Integrated to g — 23 
as expected for instance for the future BigBOSS survey, 
these counts indicat e a ~ 10% reduction compared to the 
extrapolation of the lHopkins et ail (2007) luminosity func- 
tion that was used to predict quasar nu mber counts in the 
LSST science boolfl (jAbell et all l2009f l . 



http: / / www.lsst .org/ files / docs /sciencebook 
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Fig. 13: Quasar luminosity function measurements (black circles) compared to the C09 data (blue open circles at 0.68 < 
z < 2.6). The blue dashed curve is the LF of C09 and the best-fi t model of this analysis is shown as the red curve. The 
green curves are best-fit to COSMOS data (jMasters et all 120121) at z ~ 3.2 (plain curve) and z ~ 4.0 (dashed curve). 
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Table 7: Values of the parameters of the quasar luminosity function (Eqs. [7]and|n]). The normalization of the luminosity 
function (<£>*) is given in Mpc mag -1 . The top line gives the parameters published in C09 recalculated according to 
the definition of Eq. [9] and used to initialize our fit; the bottom line is for our best-fit model on the combined data set. 



6. Conclusions 

We have designed dedicated observations to measure the 
quasar luminosity function to z = 4. The targets were se- 
lected with a technique relying on optical variability of the 
quasars, which allows both high completeness and a simple 
estimation of the incompleteness corrections. These include 
inefficiencies related to the selection technique, which we 
compute using a control sample of almost 20,000 known 
quasars, as well as instrument- related inefficiencies that we 
compute from the data. The targets were shared between 
two instruments, the Sloan telescope (though an ancillary 
SDSS-III/BOSS program) and the MMT. They yielded a 
total of 1877 quasars, divided into 436 previously known 
quasars, 1179 new quasars identified with BOSS and 262 



with the MMT. Cross-checks between the results from the 
two instruments indicate that the identification of the spec- 
tra is reliable to observed magnitudes g ~ 23. We have also 
verified that we obtained compatible number counts from 
BOSS and the MMT, in magnitude regimes where we had 
data observed with both instruments. 

These dedicated data allow us to compute reliable 
quasar counts and to derive the quasar luminosity function 
to the limiting magnitude <7dcrcd = 22.5. Relatively high 
Galactic extinction in the fields prevent us from reaching 
fainter limits. 

The quasar number counts we measure at z < 2 are 
i n agr ee ment with previous e stimates from iHopkins et all 
(|2007D or ICroom etall (|2009h . At higher redshifts, we ob- 
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Table 8: Predicted number of quasars over 15.5 < g < 25 and < z < 6 for a survey covering 10,000 deg 2 , based on our 
best-fit luminosity function. Bins are centered on the indicated magnitude and redshift values. The ranges in each bin 
are Ag = 0.5 and Az = 1. 



serve o f order 10% less quasars than assumed in lAbell et al.l 
(|2009f ). This trend is confirmed by the luminosity function 
that we fit to our data, corrected for incompleteness in the 
faintest magnitude bins using the model weighted estima- 
tor. Wh ile our best-fit mo del is in good agreement with 
that of[C room etai] ((20091) over the common range in red- 
shift 0.68 < z < 2.6, our fainter and deeper data indi- 
cate a flatter faint luminosity slope than what is predicted 
from the extrapolation of their model. This trend was al- 
re ady visible in the th e last redshift bin (2.2 < z < 2.6) 
of|C room et a l. ( 2009), and has been confirmed with the 
present work. At 3.5 < z < 4.0, our best-fit model indicates 
an excess of bright qua sars compared to the best-fit model 
of[C room et "aTT ([2009). This result needs to be confirmed 
with additional data, and would benefit from observations 
in less extincted regions of the sky in order to reach fainter 
magnitudes and increased statistics. 

Appendix 

Table [9] provides the binned luminosity function (LF) for 
the data of our analysis using the model weighted estimator 
described in Sec. [3J as plotted in Fig. Q2J We give the value 
of log<£> in 8 redshift intervals from z = 0.68 to z = 4.00, 
and for AM 9 = 0.40 magnitude bins from M g = —29 to 
M g = —21 to . We also give the number of quasars (Nq) 
contributing to the LF in the bin and the error (Alog$). 
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Table 9: Binned luminosity function using the model weighted estimator. The 1367 quasars used here are those with 
ffdered < 22.5 and resulting from the division into mutually exclusive samples as described in Sec. 14.31 and shown in 
Fig. EJJ 
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